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Summary: Dense distributions of string-like objects in material media are considered 
in terms of continuum field theory. The strings are assumed to carry a quantized abelian 
topological charge, such as the Burgers vector of dislocations in solids, or magnetic flux of 
supercurrent vortices in type-II superconductors. Within this common framework the facts 
known from dislocation theory can be extended, in appropriately modified forms, to other 
physical contexts. In particular, the concept of incompatible distortions is transplanted into 
the theory of type-II superconductivity. The compatibility law for type-II superconductors is 
derived in terms of differential forms. As a result, one obtains an inhomogeneous generalization 
of the classic Londons' equation. 
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1 Introduction 

The heuristics of the present paper can be conveniently developed by starting 
from the striking analogy between dislocations in solids and magnetic vortices 
in superconductors. While the former determine the behaviour of materials 
in plastic deformation processes, the latter are responsible for magnetic and 
transport phenomena in type-II superconductors. The physical contexts in the 
situations are indeed very different. Nevertheless, the analogy is far from be- 
ing superficial. It suffices to note common fundamental features such as the 
distinguished linear structure, quantized dynamical quantities concentrated on 
and conserved along the lines, as well as large amount of irreversibility in their 
dynamical behaviour. 

In most situations the number of vortices in a macroscopic superconducting 
specimen is very large, what enables an effective description of vortex networks 
with the aid of continuum field theories [1-3, 5-6, 8-11, 18]. 

In the present paper we shall concentrate our attention on the compatibility 
laws which govern the behaviour of the continuum fields describing the macro- 
scopic quantities. We adopt the following terminological convention. Whenever 
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more precise language concerning type II superconductivity is necessary we shall 
make use of the term fluxon in place of informal terms like flux line or magnetic 
vortex. It, however, should be contrasted with the term fluxoid, for which we 
adopt London's definition [10]. 

Thoughout the paper, a theoretical description of a string system will be 
termed macroscopic if it deals with collective quantities without paying attention 
to individual lines. The description which refers to individual strings, but not 
to individual atoms, will be called mesoscopic; the term microscopic will be 
reserved for the atomistic desription. 

2 Analogy with dislocation networks 

In general terms, the relation between magnetic flux lines and macroscopic mag- 
netic properties of superconductors is similar to the relation between dislocation 
lines and the macroscopic mechanical properties of solids. 

Burgers vector ►■ elementary flux 



dislocation lines " flux lines 



mechanical properties magnetic properties 

Later on we shall find limitations of this analogy. For the present purposes it 
turns out to be useful, and we shall try to get some profits from it. 

The continuum theory of distributed dislocations has been developed by 
Kroner in his classic work [13]. The theory has been generalized to moving 
dislocations by Kosevitch [12]. Modern geometric formulations of continuum 
theory of dislocations have been given by Kroner [14, 15], Mistura [16], and 
Trzesowski [17]. The dynamics of networks of string-like objects has been con- 
sidered recently by Rogula & Sztyren [19, 20]. 

A single dislocation line L is characterized by its Burgers vector b l L , which 
is topologically quantized and conserved along the line. On the other hand, a 
single magnetic flux line is characterized by associated magnetic flux $£, which 
also is topologically quantized and conserved along the line. An immediately 
visible difference between the two cases is the vectorial nature of b l L and the 
pseudoscalar nature of This difference, however, athough important in 

many respects, does not break the analogy. 

The presence of dislocation lines in a material continuum results in a distor- 
tion field which is well defined in the region of "good material" - outside the 
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dislocation cores. The good material is, however, multi-connected and the dis- 
locations give rise to the distortion field which cannot be represented with the 
aid of a globally defined single-valued displacement field. One of the possibili- 
ties here is to make use of multi- valued displacement fields. Another possibility 
consists in the following. Instead of a multi-valued displacement field u 1 , one 
considers a single-valued differential form du % defined by the relation 

<j> du l = 6*[C], (1) 

where the right-hand side denotes the total Burgers vector of the dislocation 
lines encircled by the contour C. The entire contour C is confined to the good 
material and is otherwise arbitrary. 

3 Macroscopic averages in string systems 

To describe effectively the material media which contain a very large amount of 
string-like objects one has the need for a fully macroscopic theory which operates 
with quantities describing large collections of such objects without paying much 
attention to an individual string. 

3.1 String density and string current 

The placement and movement of an individual line segment at the mesoscopic 
level can be conveniently described with the aid of the following quantities, valid 
generally for networks of moving lines L = L{t) 

a 3 (x | L) = J 5(x-x')dx' j (2) 

and 

L k (x,t | L(t)) = / <5(x - ^)e k:j iv l (^)dx\ (3) 

JL(t) 

where v (x ) denotes the velocity of the line element at x . Only the components 
of velocity perpendicular to the line element contribute to the string current (3). 

Note that the above defined quantities, which may be interpreted as the ori- 
ented string density and oriented string current, depend solely on the orientation 
and geometry of the lines, and are independent of any physical characteristics 
of the strings. For the string density and string current we obtain the following 
balance equations 

^ + ^V,=0, a* i< = 0. (4) 
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It is expedient to notice that, in the case of point-like particles, the analogue 
of the first of the above equations expresses the conservation of the number of 
particles. 

For networks of dislocations in crystals and magnetic flux lines in supercon- 
ductors, the related quatities are defined in the following subsections. 

3.2 Dislocation networks 

Notation: b l L represents the Burgers vector associated with the line segment L. 

a ij (x | L) = b l L a j (x \ L), (5) 

< es >) = £^'(*|£)> (6) 

L 

= (7) 
J\ meso (x) = J2 b W(x,t\L(t)), (8) 

L(t) 

J i k (x,t) = (J i k mes °(x,t)) av . (9) 
The dislocation density a y and the dislocation current J l k satisfy the equations 

!^' + ^A = o, <i = o (io) 

which express the Burgers vector conservation laws. 

3.3 Magnetic flux networks 

Notation: and <&o represent the quantized flux associated with the line 
segment L and the elementary flux quantum, respectively. Typically <&l = ±$o- 

Q j (x | L) = § L a j {x | L). (11) 

& meso (x) = J2 Q3 ( x \ L )> ( 12 ) 

L 

= (OLsoWU, (is) 
Kr s °(x) = Y,*LH(x,t\L(t)), (14) 

L(t) 

Ki(x,t) = (K™ s °(x,t)) av . (15) 
The density 6^ and the current Ki satisfy the conservation equations 

*e< + e«^ = 0, 6^ = (16) 
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which are the analogues of eqns. (10). The quantity Q l represents the contri- 
bution of the magnetic flux lines (lluxons) to the macroscopic induction field 
B % \ therefore it will be called the (oriented) fluxon density. The quantity Kj 
represents the macroscopic fluxon current density resulting from the motion of 
fluxons. On the other hand, the flux current produces some electric field, cou- 
pled to the moving magnetic flux through Maxwell equations. In consequence, 
we arrive at the following interpretation 

B str = 9, E str = -K, (17) 
c 

where B str and E str stand for the contribution of strings to the electromagnetic 
fields. 

Several approaches to the continuum description of superconductors contain- 
ing dense distributions of vortex lines are known from the literature; see e.g. 
London & London [9, 10], Laue [8], Zhou [11], Abrikosov [1, 2], Anthony & 
Seeger [3], Chapman et al. [5, 6], Rogula [18]. Some of them can be derived 
with the aid of particular assumptions concerning the macroscopic averages. 
Note that, in general, due to statistical correlations between mesoscopic string 
configurations and their motions, the inequality 

(a j v l ) ^ (a j )(v l ) (18) 

holds. There are, however, important special cases, such as coherent fluxon flow, 
or else uncorrelated randomness of a and v, when both sides of the formula 
(18) are identical to a good degree of approximation, and the corresponding 
assumption of the vortex-density model [6] is justified. The fluxon current K 
and the fluxon density O are then related through the macroscopic average of 
the fluxon's velocity by the equation 

K + vAQ = 0, (19) 

valid in this special case. 

4 The compatibility equations 

4.1 The distortion field and dislocation density 

The presence of continuously distributed dislocation lines in a material con- 
tinuum causes an incompatible macroscopic distortion of the medium. The 
situation may be briefly sketched as follows. While the multi-valued mesoscopic 
displacement fields, mentioned in Section 2, are rather hard to be averaged, and 
would lead to ill-defined macroscopic diplacement fields, the single- valued differ- 
ential forms allow for unequivocal averages. Therefore, to describe incompati- 
ble distortion fields at the macroscopic level, instead of ill-defined displacement 
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field u % one considers the well-defined differential form du % undestood as the 
macroscopic average of the mesoscopic differential form (1). The corresponding 
continnum field equations, expressed in terms of the material distortion /3* and 
the dislocation density a 11 , take the form 

e lk > /3i k =a il , a u t=0, (20) 

and can be obtained from the expression 

du l (x,t) = /3 i j (x,t)dx j (21) 

by taking into account the equations (1) and (7). Keeping that in mind, we will 
search for an analogue of the above relations in the framework of superconduc- 
tivity. 



4.2 The fluxoid density and superconducting vortices 

Let us consider the expression 

h = - M) -^- c \4>\ 2 A k . (22) 

for the supercurrent density (for simplicity we start from an isotropic super- 
conductor) with the following notation: ip — ip(x,t) represents the Ginzburg- 
Landau order parameter, -0 its complex conjugate, A k = A k (x,t) the electomag- 
netic vector potential, e* and m* denote the electric charge and the effective 
mass of a Cooper pair, respectively. We substitute 

V> = I ip I z, (23) 

where | z \ = 1 so that z = z(x, t) equals the phase factor of the order parameter 
field. The eqn. (22) may then be rewritten as 

3k = — I ^ I 2 {^(zz, k -z, k z) A k ). (24) 

Now, let us introduce the following differential form: 

d(j) = f —izdz = i dzz, (25) 
which represents the phase differential. With the aid of this form we can write 

j k dx k = | V | 2 (hdcf> + -A k dx k ) (26) 

m* c 

or, taking into account that in the superconducting region outside the vortices 
4> + 0, 

-^ =(4+ S^f^- (27) 
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The right-hand side of the above equation represents the fluxoid density. After 
integration over an appropriately smooth surface S, we obtain 



he 

e* 



!0 = *[C], (28) 
where $[C] equals the total fluxoid encircled by the contour C = dS, 



$[C] = f e'dflj. (29) 
Js 



In consequence, passing to the differential relation, we obtain the analogue 
ofcqn. (20) 

^ + $|f]2)>i = 0i > ^ = (30) 

or equivalently 

B i +e ilk { ^ jkU = Q i t (31) 



with n s = \ ip | 2 . The left-hand side of this equation corresponds to the gauge 
invariant Londons' equation 

B i +e ilk { J^ jk) (32) 

with the charge e and mass to replaced by the Cooper pair charge e* and 
effective mass to*, respectively, and with | V 1 2 interpreted as n s . Hence one 
can see that the incompatibility of the superconducting order parameter due to 
continuously distributed flux lines modifies Londons' relation (32) between the 
magnetic induction field and the supercurrent density. 

Due to the electromagnetic character of the quantities involved, an electric 
counterpart of eqn. (31) has also to be valid. In fact, taking into account the 
flux flow balance (16) and making use of the Maxwell equation 



Jjk 

we obtain 



1 dB l 



' J " Ek ' j + ~ c ~dt = °' (33) 



e ilm (^-A^-j m ) + -K m -E m ) tl = 0. (34) 
ot e* z n s c 

The above equation can be conveniently integrated to the form 



dm* 1 
Ek = - X , i + m {^- ji ) + -K i . (35) 



Note, however, that the integrated relation (35) becomes algebraically indepen- 
dent of the equations derived up to this point. This is due to the new gradient 
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term x,i which depends on specific features of the configuration under consid- 
eration. 

The equations given so far in the present section are valid literally for 
isotropic superconductors. The generalization to anisotropic superconductors 
can be, however, performed in a straightforward manner by applying the rou- 
tine procedure based on substitution of the effective mass tensor m*j in place 
of the scalar mass m* . The tensor m*j is real and symmetric by definition, and 
positive-definite by assumption; in particular, it has a well defined inverse which 
can be substituted to eqn. (22). As a result, the equations (31) take the form 



47T 



and 



The tensor 



B l + —e il i(\%j k ), l = e i 



_ 47T 8 ,,o .u, 1 T . 



A 2 = <? m *jk 

° k Ait e* 2 n s 



(36) 

(37) 
(38) 



generalizes the (squared) Londons' penetration depth \\. 

5 Conclusion 

The heuristic formulation of the analogy between dislocations in solids and flux 
lines in superconductors can now be stated in a more precise way: the diagram 
given in the Introduction should be complemented by the following scheme: 



distortion field 



dislocation density 



dislocation current 



fluxoid density 



fluxon density 



fluxon current. 



Due to the presence of extra quantities - the fluxon density and the fluxon 
current - the above given macroscopic equations, even when augmented with 
material constitutive relations, are incomplete. They can be completed by a va- 
riety of mathematical models stated in terms of the fields <d l and Kj and taking 
into account the string kinetics. In this way one can, for instance, reproduce 
the behaviour of superconductors described by critical state models, such as 
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Bean's [4] or Kim & Stephen's [7] ones. On the other hand, by making use of 
the above stated analogy, one can adapt a selection of models conceived in the 
theory of dislocations in order to describe the irreversible behaviour of dense 
string systems. 
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